The Fine Structure Lines of Hydrogen in HII Regions 
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The 2si/2 state of hydrogen is metastable and overpopulated in HII regions. In 



addition, the 2p states may be pumped by ambient Lyman-a radiation. Fine structure 
transitions between these states may be observable in HII regions at 1.1 GHz (2si/2 — 
^ I 2pi/2) and/or 9.9 GHz (2si/2 — 2^3/2), although the details of absorption versus emission 

are determined by the relative populations of the 2s and 2p states. The n = 2 level 
populations are solved with a parameterization that allows for Lyman-a pumping of 
j_j ■ the 2p states. The Lyman-a pumping rate has long been considered uncertain as it 



involves solution of the difficult Lyman-a transfer problem. The density of Lyman-a 
photons is set by their creation rate, easily determined from the recombination rate, 
and their removal rate. Here we suggest that the dominant removal mechanism of 
Lyman-a radiation in HII regions is absorption by dust. This circumvents the need to 
solve the Lyman-a transfer problem, and provides an upper limit to the rate at which 
the 2p states are populated by Lyman-a photons. In virtually all cases of interest, 
the 2p states are predominantly populated by recombination, rather than Lyman-a 
pumping. We then solve the radiative transfer problem for the fine structure lines in 
the presence of free-free radiation. In the likely absence of Lyman-a pumping, the 
251/2 ~^ 2pi/2 lines will appear in stimulated emission and the 2si/2 2p^/2 lines in 
absorption. Because the final 2p states are short lived these lines are dominated by 
intrinsic linewidth (99.8 MHz). In addition, each fine structure line is a multiplet of 
three blended hyperfine transitions. Searching for the 9.9 GHz lines in high emission 
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measure HII regions offers the best prospects for detection. The Unes are predicted 
to be weak; in the best cases, Hne-to continuum ratios of several tenths of a percent 
might be expected with line strengths of tens to a hundred mK with the Green Bank 
Telescope. Predicted line strengths, at both 1.1 and 9.9 GHz, are given for a number 
of HII regions, high emission measure components, and planetary nebulae, based upon 
somewhat uncertain emission measures, sizes, and structures. The extraordinary width 
of these lines and their blended structure will complicate detection. 

Subject headings: ISM: HII regions — ISM: lines and bands — line: formation — 
planetary nebulae: general — radio lines: ISM 

1. Introduction 

Transitions between fine-structure sublevels in atomic hydrogen have never been detected as- 
tronomically. With advances in radio astronomy, however, it now appears that transitions between 
the 2si/2 and 2pi/2, and 2si/2 and 2p3/2 levels may soon be observable in HII regions. Figure 1 
shows the fine and hyperfine structure of the n = 2 level of atomic hydrogen. The fine structure 
splitting into the 2si/2, 2pi/2i and 2p^/2 levels is caused by a combination of spin orbit coupling, 
relativistic effects and the Lamb shift. Significantly, the 2si^2 state is mctastable owing to the 
transition rules for angular momentum, i.e. the 2s-ls transition is forbidden. Hydrogen atoms in 
the interstellar medium are removed from the 2si/2 state principally through 2-photon emission to 
the ground state (Breit and Teller 1940) with decay rate per atom, ^27 = 8.227 sec~^ (Spitzer and 
Greenstein 1951), resulting in optical continuum radiation. Because this rate is some eight orders of 
magnitude slower than Lyman decay, it may seem reasonable to expect that the 2si/2 state will be 
overpopulated relative to other excited states under photoionization equilibrium. It must be noted, 
however, that trapped Lyman-a radiation in an HII region will pump the 2p states. In addition, in 
dense HII regions collisions, primarily with ions, will transfer hydrogen atoms between the 2s and 
2p states. These effects must be taken into account in order to determine the expected strengths 
and signs (absorption versus stimulated emission) of the 2si/2-2pi/2 (1-1 GHz) and 2si/2-2p3/2 (9-9 
GHz) lines. 

Apparently, Wild (1952) was one of the first to investigate the astrophysical implications of 
these lines, suggesting that the 2si/2-2p^/2 liiie might be observed in the sun. Purcell (1952), how- 
ever, noted that collisions in the dense solar atmosphere would equilibrate the 2s and 2p populations 
making detection unlikely. Townes (1957) suggested that an overpopulation of the metastable 2s 
states may lead to detectable lines in the interstellar medium. Shklovski (1960) raised the possi- 
bility of Lyman-a pumping of the 2p states but discounted its sufficiency for producing observable 
2p 2s transitions. Pottasch (1960) argued that the Lyman-a depth in HII regions is essentially 
infinite such that virtually all downward Lyman-a transitions are balanced by a Lyman-a absorp- 
tion, resulting in an overpopulation of 2p states relative to 2s at densities below about 10^ cm~^. 
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Fig. 1. — Fine and hyperfine structure of the n = 2 level of hydrogen. Allowed hyperfine transitions 
are shown as solid vertical lines, with arrows denoting predominant transition directions if the 
metastable 2si/2 states are overpopulated with respect to the 2p states. (See text.) Hyperfine 
transition frequencies are given in MHz. The quantum number F corresponds to total angular 
momentum, i.e., orbit + electron spin + nuclear spin. Note that level separations are not shown 
to scale. 
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Field k. Partridge (1961) further considered the case in which Lyman-a radiation is effectively 
destroyed by collisional conversion of 2p states to 2s (followed by 2-photon decay to the ground 
state). The consequent overpopulation of 2p states then led to the prediction that the 9.9 GHz 
lines would appear in stimulated emission in Hll regions. (This also implied that the 1.1 GHz lines 
would appear in absorption.) Myers k, Barrett (1972) used the Haystack 37-m radiotelescope with 
a 16-channel filter bank spanning 1 GHz in frequency to search for the 9.9 GHz lines, but did not 
find any lines in excess of 0.1 K antenna temperature. Ershov (1987) estimated the strengths of 
both sets of lines in Orion A and W3(0H) under the assumption that the 2p states are negligibly 
populated with respect to 2sxi2- 

The major variation in the predictions arises from the uncertain density of Lyman-a radiation 
in HII regions, which in turn is determined by the processes that destroy Lyman-a radiation. In 
addition to the conversion of 2p to 2s states via collisions, other processes are likely to be important. 
Resonantly scattered Lyman-a photons undergo a diffusion in frequency and may thus acquire a 
significant probability of escaping the HH region due the much lower optical depth in the line wings 
(Cox & Matthews 1969; Spitzer 1978). Systematic gas motions within the region (e.g. expansion) 
can also contribute to the migration of Lyman-a photons away from line center. The dominant 
removal mechanism, however, is absorption by dust within the HII region (Kaplan h, Pikelner 1970; 
Spitzer 1978). This last fact greatly simplifies the estimation of the Lyman-a density. 

In this paper, we obtain formulae for the 2s and 2p populations in an HII region under a 
parameterization that characterizes the rate of 2p production by Lyman-a absorption in terms 
of the rate of 2p production through recombination (Section 2). Existing models for dust in HII 
regions and planetary nebulae are then used to place tight constraints on the Lyman-a density 
and thus the 2p production rate by Lyman-a photons (Section 3). The radiative transfer problem 
for the fine structure lines is then solved for a uniform region (Section 4). We obtain approximate 
predicted line strengths for various HII regions, compact components, and planetary nebulae under 
the justified assumption that Lyman-a pumping of the 2p states is negligible (Section 5). High 
emission measure HII regions and/or planetary nebulae are most likely to yield detectable lines. 



2. Population of the 2s and 2p States 

Spitzer h Greenstein (1951) estimated that between 0.30 and 0.35 of all recombinations in an 
HII region reach the 2s state. Including the effects of collisions and the 2-photon decay rate, the 
rate equation for the 2s state is 

fauiUe + Cpsnin2p = ^27?^2s + Cspnin2s, (1) 

where / is the fraction of recombinations producing the 2s state and a is the recombination co- 
efficient. The collision coefficients, Cpg and Cgp, give the appropriate rates at which atoms are 
transferred from 2p to 2s and from 2s to 2p, respectively. These rates are dominated by collisions 
with protons, but including electron collisions (and assuming that Ug = ni and T = 10^ K), we have 
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Csp = 5.31 X 10 ^ cm^ s ^ (Seaton 1955). Since kT is much greater than the energy separations 
between the 2s and 2p states, the reverse rate is determined by just the ratio of statistical weights, 

i-e- ^ 

Cps = Cgp = '^Cgp- (2) 
92p 3 

Similar considerations apply to the 2p states. The fraction of recombinations reaching 2p is 
(1 — /). [It is worth noting that processes such as Lyman-/? absorption, which ultimately produce 

the 2s state through Is ^ 3jo — > 2s, have already been taken into account in the calculation 
of / (Spitzer &: Greenstein 1951).] Radiative decay occurs through both spontaneous Lyman-a 
emission (with A21 = 6.25 x 10^ sec^^), as well as stimulated emission. We must also include 
radiative excitation from the ground state via Lyman-a absorption. Including collisions which 
couple to the 2s states, the rate equation for 2p then is 

{1 — f)anine + Cspnin2s + c J n^J^^dv J Bi2{v — v')nu'dv' 

= Min2p + Cpsnin2p + c J nf^^dv J B2i{i' - u')ny'dv' , (3) 

where n^^"^ and rii^^^ are the densities of atoms in the Is and 2p states per radial velocity interval, 
measured in frequency units; and ni,i is the photon density per frequency interval (in cm~^ Hz^^). 
Now, 



B2i{u-u') = ^A2iL{v-v') (4) 
and Bi2{v-u') = W2i{v-v'), (5) 
where is the Lyman-a frequency, and L{u — v') is the Lorentz line profile. For a thermal gas, 

4^^) = -^^e'^, (6) 

where (Ai^)d is the Doppler width. For T = 10^ K, {Av)d = 1-29 x 10^^ Hz. Clearly thermal 
widths will completely dominate the natural (Lorentz) line width of Lyman-a, and thus we can 
replace L{i' — v') by the Dirac delta function. We then obtain for the absorption term: 

Of course, a similar calculation can be carried out for the stimulated emission term. 

Because the nebula is optically thick in Lyman-a, must be considered carefully. Were 
escape from the nebula the primary removal mechanism, then a steady state would result in which 
photons are created near line center, diffuse in frequency through resonant scattering and are 
effectively removed far out in the wings. The decrease in both creation rate and diffusion rate 
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with frequency offset results in a nearly flat distribution (n^) out to the approximate frequency 
at which photons freely escape, beyond which it drops sharply (Capriotti 1966). We write this 
frequency offset as w{Aiy)D, where it; is a dimensionlcss parameter expressing the offset in terms of 
the Doppler width. We shall therefore assume that rii^ = nLya/['2iv{Ai')£)] for [v — i^dI < w{Ai')£), 
and n^, = for — > w{Aiy)£), where niya is the density of Lyman-a photons. The 2p rate 
equation then becomes 

n ^^ , 3c^ A21 erf(w) 

(1 - JjaniUe + Cspn2sni + — — oTTT^ riLyanis = ^21^2^. (8) 

Spontaneous emission completely dominates the depopulation of the 2p state and thus the other 
two terms that appeared on the right hand side of equation (3) have been dropped. Stimulated 
emission is negligible in comparison with spontaneous emission for any reasonable value of riLya. 
Indeed, equality of spontaneous and stimulated emission would imply a radiation pressure (due to 
Lyman-a) many orders of magnitude in excess of the thermal gas pressure. Also, the 2p 2s 
collision rate is negligible for any reasonable value of rij. [We include these collisions terms insofar 
as they populate the 2s state, however. Sec equation (1).] 

As we shall see, it is quite plausible that Lyman-a photons are absorbed by dust before 
significant frequency diffusion occurs. In this case n,y will simply reflect the thermal distribution of 
atoms. The resulting 2p rate equation is identical to that given above provided we replace eif{w)/w 
with 2sj2fK. 

At low densities, Lyman-a photons are created at the rate at which recombinations lead to 2p 
states, i.e. {1 — f)anine- At high densities (n^ > 10^ cm~^) 2s states may be collisionally converted 
to 2p states leading to additional Lyman-a photons, and thus the Lyman-a creation rate could be 
as large as anj^e, the total recombination rate. We therefore define the Lyman-a lifetime as 

thya — ) (9) 

rariiUe 

where 1 — / < r < 1. Since 2p states decay to the ground state much faster than they could be 
collisionally converted to 2s, all recombinations to 2p are regarded as producing a Lyman-a photon. 
Thus, r could never be smaller than 1 — /. 

We define 

3c^ A21 r erf(w) 
S = — 2 /A \ n 7\ XnntLya, (10) 

le-TTI/'f (Az/)i3 (1 — /) W 

where x = nis/nff, and nn is the total hydrogen density (atomic plus ionized). The rate equations 
[(1) and (8)] can then be solved for the 2s and 2p populations: 



a nine[fA2i + (1 - /)(1 + 5)Cp. n,] 
^21(^27 + CgpUi) — CspCpsuf 



= ^— n , ^ ^ ^ ^ o (J-lj 



_ an,ne[(l - /)(1 + S)A2^ + {I + S - /5)C,pn,] 

A2l{A2j + CspUi) - CspCpsUf 
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The dimensionless parameter S gives the rate at which 2p states are produced by captured Lyman-a 
radiation in terms of the 2p creation rate from recombination (at low densities). 

The ratio n2p/ (3n2s) provides a determination of the relative importance of Lyman-a pumping 
of the 2p states and whether a fine structure line is expected to appear in absorption or (stimulated) 
emission. If n2p/{3n2s) > 1, then the 9.9 GHz, 2si/2-2p3/2 line will appear in emission, and the 1.1 
GHz, 2si/2-2pi/2 line will appear in absorption. Of course, n2p/(3ra2s) < 1 implies the opposite. 
(This assumes that the two 2p states are populated according to their relative statistical weights. 
To evaluate cases in which n2p/(3n2s) is of order unity, it would be necessary to separately account 
for the rates at which 2pi/2 and 2^3/2 states are created and destroyed, including the coUisional 
rates coupling these states.) Prom equations (11) and (12) it follows that n2p/ (3^125) > 1 if > Scrit 
where, 

5c.^t = -r^4^ = 1.14xl0«. (13) 

(The numerical value of Scrit corresponds to / = 1/3.) Because the 2p states naturally decay 

about 10^ times faster than the 2s states, population equality thus requires a pumping rate some 8 
orders of magnitude faster than the approximate rate at which 2s and 2p states are formed through 
recombination. 



3. The Lyman-a Density in HII Regions 

Determination of the Lyman-a density in HII regions is a complicated transfer problem. It 
is likely, however, that the dominant mechanism for removing Lyman-a photons is quite straight- 
forward, i.e. absorption by dust (Kaplan Sz Pikelner 1970; Spitzer 1978). Thus, we eschew the 
noncoherent radiative transfer problem and find the upper limits to t^ya and S set by absorption. 
Other competing removal processes would reduce the lifetime, and therefore density, of Lyman-a 
photons, resulting in a lower value for S. 

Using a silicate-graphite model for dust in HII regions (Aannestad 1989), with a dust-to-gas 
ratio of 0.009, the extinction at Lyman-a can be shown to be Nh/{5A x lO^^'cm^^) mag, where Nh 
is the column density of hydrogen. The albedo for this mixture is about 0.4 at Lyman-a (Draine 
Sz Lee 1984). The lifetime of Lyman-a photons against absorption by dust can then be calculated 
as tLya = (3.3 X 10^° cm~^ s)/nH- We then find that 

S = 4.2X10^X ,/,,"^'"" . (14) 
(1 - /) w 

Since x << 1 throughout most of the volume of an HII region (Osterbrock 1989), we conclude that 
S « Scrit and that the 2s state is overpopulated relative to the 2p states. 

In the harsh environments of planetary nebulae dust might be destroyed by shocks or hard 
UV radiation, or possibly separated from the ionized gas by radiation pressure (Natta Sz Panagia 
1981; Pottasch 1987). Abundance measurements indicate, however, that various heavy elements are 
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depleted from the gas phase in the ionized regions of NGC 7027 (Kingdon &: Ferland 1997) and NGC 
6445 (van Hoof et al. 2000), suggesting that dust has not been destroyed in significant quantities. 
Additionally, planetary nebulae frequently exhibit a mid-IR spectral component characteristic of 
warm dust heated by the intense radiation field within the ionized region (Kwok 1980; Hoare 1990; 
Hoare et al. 1992). Summarizing these results, Barlow (1993) has argued that dust is a common 
constituent of the ionized zones of planetary nebulae, albeit with dust-to-gas ratios about an order 
of magnitude or more below that of the general ISM. 

Middlemass (1990) has modeled NGC7027 and finds an extinction optical depth of about 0.17 
at 500.7 nm. For a uniform model the column density of hydrogen in the ionized zone is about 
3.5 X 10^^ cm~^ (Thomasson & Davies 1970), along the radius of the nebula. For the graphite 
dust model used by Middlemass (1990) the extinction at Lyman-a then is Nh/{1.?) x lO^^cm"^) 
mag; about a factor of 20 smaller than that in a general HII region described above, and roughly 
consistent with the smaller dust-to-gas ratio (« 7 x 10~^) in this object (Barlow 1993). For an 
albedo of 0.4, the Lyman-a lifetime is tLya = (7.6 x lO^^cm"^ s)/nH, and thus 

S« 9.8X10^X7^^. (15) 
(1-/) w 

If we assume that absorption by dust is the dominant process limiting the Lyman-a density, make 
the replacement eif{w)/w 2^j2fK^ and let r f« 1, we obtain S = 2.3 x lO^X- Since it is quite 
unlikely that x is as large as 0.05, as required for S ~ Scrit^ we conclude also in this case that the 
2p states arc underpopulated relative to the 2s states. 

Of course, the above estimates are upper limits to the Lyman-a lifetime, as other processes 
(described in Section 1) may contribute to the removal of these photons. Thus, values of S estimated 
in this way are upper limits. In addition, if escape is important then the relevant value of eYi{w)/w 
would be smaller than that substituted above. 

It should also be noted that Lyman-a radiation contributes to heating the dust in the ionized 
region. The energy absorbed is then reradiated in the IR. In the case of NGC 7027 a mid-IR spectral 
component is evidently due to dust with temperature T^i ~ 230 K comixed with the ionized gas 
(Kwok 1980). If the heating is dominated by Lyman-a radiation, then 

riLyahvLcQabsiLya) = 4 {Q{a, T^)) aT^, (16) 

where Qabs{Lya) is the absorption efficiency at Lyman-a and {Q{a,Td)) is the Planck-averaged 
emissivity (Draine & Lee 1984), a function of grain size a and temperature T. Large grains, being 
efficient emitters in the IR, will yield a larger estimate of niya- Thus, we assume 1 ^um graphite 
grains for which Qabs{Lya) « 1 and (Q(l ^um, 230 K)) f« 0.05 (Draine & Lee 1984). Then, 

n,„,|^ 'Q'"--^^'"q> ,6.5xl0^c.n-». (17) 

hVLC Qabs 



We then find 



3c^ A21 1 erijw) x ^Lya _ ip erf jw) x 

167rul{Ai.)Da{l-f) w (1 - x) " w {1 - xY ^ ' 
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where we have taken rij = (1 — x)''T'H (since the fraction of atoms in excited, bound states is 
negligible). To obtain the numerical value given above we used the recombination coefficient in 
the density bounded case with T = 10^ K (Osterbrock 1989). Assuming Lyman-a radiation is 
primarily removed through absorption by dust (in which case erf(it;)/u) is replaced by 2y^2/7r), 
then S < Scrit unless x exceeds 2.4 x 10~^, which is unlikely (Osterbrock 1989). It should also be 
noted that other sources, such as continuum radiation, are likely important in heating the dust, 
thereby reducing further the implied value of ULya (and therefore S) . Finally, the 1 /im grain size 
assumed here is probably an overestimate implying that niya is also overestimated. 



4. Radiative Transfer of the Fine Structure Lines 

Evidently, the 2s state is overpopulated relative to 2p, thus the fine structure transitions will 
proceed from 2sx/2 to 2p3/2 via absorption and to 2pi/2 via stimulated emission. Although the 
2p populations are probably negligible, they will be included in the radiative transfer calculation. 
The distribution of 2p states between 2pi/2 and 2^3/2 may deviate somewhat from the statistical 
weights (1/3 and 2/3, respectively), in part because the separate collisional rates from 2s are not 
proportional to the statistical weights. Since the 2p population is most likely negligible, a detailed 
calculation of the distribution between 2pi/2 and 2^3/2 states will not be carried out here. Rather, 
the fractional populations of 2px/2 and 2^3/2 will be parameterized as /3a/3 and 2/36/3, respectively. 
If /3a = /Jft = 1, then these states are populated according to their statistical weights. There is the 
obvious constraint that (3a/S + 2/3b/3 = 1. 

The fine structure transitions are allowed electric dipole transitions and the corresponding rates 
may be computed in a straightforward manner (Bethe &; Salpeter 1957), giving Aa = 1.597 x 10"^ 
sec"-*^ (2si/2-2pi/2) and Af, = 8.78 x 10""^ sec""-*^ (2p3/2-2si/2)- The absorption coefficient (valid for 
either transition) is 

8ttu^ g2s 3 

where g is the degeneracy of the final state (2 for 2pi/2j 4 for 2^3/2; g2s = 2). The — sign corresponds 
to transitions to 2pi/2', the + sign to transitions to 2^3/2, and P is either /?„ or /^t, respectively. 
Either final state quickly decays to the ground state via Lyman-a with rate A21. The natural line 
width of the 2s-2p transitions is therefore dominated by the rapid decay of the 2p state and is 
r = A2i/2'K = 99.8 MHz. For a Lorentzian profile, 

_ A{r/2.) 
'^"(---/)^ + (W ^ ^ 

where A is either Aa or A}, and Uf is the frequency of the fine structure transition. Prom equations 
(11) and (12) we find 

^2 - -n2 - fo^^i^e I3S ^ ^^^^ 

3 (^27 ~l" CgpTli) Scrit 



'^^ = =^ Q^, 2 ... Ai,{n2s - -^ri2p), (19) 
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where we have assumed that the Lyman decay rate ^421 always dominates the colhsion rate CspUi 
The absorption coefficient at Une center then is 



g 2 A fanjUe ,^ _ _PS_ 

Scrit 

This will be scaled according to the free-free absorption coefficient (Altenhoff et al. 1960), 



ff 0-212n,ne 
'^'^ ~ ^2.12^1.35 ■ ^'^'^) 

The ratio of the line optical depth (at line center) , r^"^ , to the free-free continuum optical depth, 
ri^ , then is 

max jy- f3q 

where K = —3.0 x 10~^ sec~^ for the 2si/2 transition and K = 0.41 sec~^ for the 

2si/2 — 2^3/2 transition. We assumed T = 10^ K, and that all ionizing photons are captured by 
the HIT region (Osterbrock 1989). 

The equation of transfer is 

^ = {-KfJ -K,)h + JiL (25) 

Note that the line produces absorption or stimulated emission only (through Ky), but does not 

produce significant spontaneous emission. Thus, the emissivity is completely dominated by the 
f f 

free-free emissivity, Jo . For an HII region uniform along a ray path the solution is 

= ^ 1 + R • 

where the result has been expressed in terms of brightness temperature Tj,. The line strength (in 
K) divided by the continuum brightness temperature is 



AT. 1 - 



Tb,cont 1+R e^^^ - 1 
In the optically thin limit {tI^ << 1) expanding (to second order) gives 

ATb 1„ „ 1 



-1. (27) 



Rrlf = --C-. (28) 



Tb^ccmt 2 2 

This result is in agreement with that obtained by Ershov (1987) for the optically thin limit. In the 
optically thick limit (r/-^ >> 1), 

ATb -R , , 
— = . (29) 

,cont 

1 + R ' ' 

The lines do not disappear in the optically thick limit as long as the 2s and 2p states are not in local 
thermodynamic equilibrium. (See also Ershov 1987.) The 10^ K microwave radiation field is not 
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sufficient to overcome tfie rates described in Section 2 and tfius it can not establish equilibrium. For 

example, in the likely case where n2p/{2>n2s) « 1, the 2p levels are drained by Lyman-a emission 
faster than the microwave radiation field can return these states to 2s via either absorption or 
stimulated emission. (See also Section 2.) 

Each fine structure level is split into two hyperfine levels as shown in Figure 1. The allowed 
transitions are also indicated. Both fine structure lines are split into three hyperfine components. 
The relative intensities of the components can be calculated from the appropriate sum rules (So- 
belman 1992). For the 910, 1088, and 1147 MHz lines the ratios are 1:2:1. The 9852, 9876, and 
10030 MHz lines appear in the ratios 1:5:2. 

5. Prospects for Detection 

The 9.9 GHz transitions are intrinsically about three orders of magnitude stronger than the 
1.1 GHz transitions. The solution to the radiative transfer equation indicates that, in general, the 
line brightness temperature ATf, grows as the square of the free-free optical depth in the optically 
thin limit (equation 28). One factor comes from the growth of free- free emission, which effectively 
forms the "background"', and the other from the proportionate growth of the line absorption. In 
the optically thick limit, the growth in line temperature saturates (equation 29). Most HH regions 
of interest are optically thin at 9.9 GHz and optically thick at 1.1 GHz. In general, the higher 
optical depth at 1.1 GHz docs not significantly offset the relative weakness of these transitions. We 
find, therefore, that the 9.9 GHz lines offer the better prospects for detection. 

Table 1 gives estimates of the line strengths for various HII regions and components, and 
planetary nebulae. It was assumed that S* = 0, in which case the 9.9 GHz lines would appear in 
absorption and the 1.1 GHz lines in stimulated emission. Most of the entries in Table 1 correspond 
to high emission measure components in HII regions, which are listed according to the nomenclature 
of the original references. The published emission measure values E were used to calculate the free- 
free optical depth and the continuum brightness temperature (assuming in most cases Tg ~ 10^ K). 
The line-to-continuum optical depth ratios R arc calculated from the published electron densities 
(equation 24). Notably, our predictions for both the 1.1 GHz and 9.9 GHz lines from Orion A 
(M42) are in good agreement with those of Ershov (1987). 

The estimates in Table 1 take into account the distribution of the line strength over three 
hyperfine components making up each fine structure line and the consequent line blending. At 
9.9 GHz, the strongest line (9876 MHz) will be blended with the weakest fine (9853 MHz). The 
peak temperature occurs at 9874 MHz and is 75% of that calculated using equation 26 for a 
single (fictitious) fine structure line. The 10030 MHz line will be somewhat distinct with a peak 
temperature equal to 32% of that predicted from equation 26, including contributions from the 
wings of the 9852 and 9876 MHz lines. The situation at 1.1 GHz is similar. The 1088 MHz hue 
will appear blended with the 1147 MHz line, with a peak value of 63% at 1093 MHz. The 910 
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MHz line will remain distinct with a peak value of 30%, including line wing contributions from the 
other two lines. The estimates in Table 1 give the line temperature at the peak of the brightest 
spectral feature in either the 9.9 GHz or 1.1 GHz blended multiplet. At both frequencies, the 
resulting profile of a strong blended feature and a distinct weaker line, combined with Lorentzian 
line shapes, will provide a unique detection signature. Figure 2 shows the predicted appearance of 
the blended multiplet around 9.9 GHz. 

The Green Bank Telescope (GET) provides a realistic hope of detection of the 9.9 GHz lines. 
Having a clear aperture it should be relatively free of standing waves in the antenna structure. 
Thus, it should be possible to search for very broad spectral features over bandpasses as large as 
800 MHz. Nevertheless, the high emission measure HIT regions that are likely to show the lines 
are typically quite compact. Therefore, the estimated line antenna temperatures (ATa)9.9 include 
the effects of dilution in the l'.21 GET beam. The estimated values of AT^ at 1.1 GHz were 
calculated using either the 11' GET beam (subscript G), or where appropriate, the Arecibo 4'. 3 
beam (subscript A). 

Quite possibly the observational sensitivity will be limited by systematic effects in which the 
continuum antenna temperature is modulated by frequency dependent gain variations which are not 
entirely removed by calibration procedures. Thus the limiting factor may be the line-to-continuum 
ratio, AT/T, which is also tabulated for both lines. 

The estimates given in Table 1 are only approximate. The underlying observations are biased 
in favor of component sizes to which various interferometer arrays are sensitive. In complex HH 
regions structures are present on a range of scales, down to very high emission measure arc-second 
scale components (Turner &; Matthews 1984). The components included in the table (typically a 
few tens of arc seconds in size) were selected because they contribute significantly to the total flux 
in the GET beam. The more compact, arc-second scale components typically yield larger values 
of AT/T (despite densities > 10^ cm~^ and consequent collisional de-excitation), yet they tend to 
contribute relatively little to the total flux in the GET beam. Conversely, extended, low emission 
measure components in the beam will add to the antenna temperature while contributing little line 
absorption. It should also be noted that the presence of structures having a wide range of densities 
implies that the uniform model considered in Section 4 is an oversimplification. Inhomogeneous 
structure, including clumping, in the emission regions would imply that the emission measure 
estimates in Table 1 represent an average over the surface of the source. In the optically thin case, 
i.e. most sources at 9.9 GHz, a redistribution of emission measure, and therefore optical depth, 
will tend to strengthen the overall estimated line strength AT5 due to its r^-dependence (equation 
28). Collisional de-excitation in denser regions of the source, however, will tend to reduce the line 
optical depth (equation 24). Eecause R is density dependent, particularly for densities above about 
1.5 X 10^ cm~^, the brightness temperature for a homogeneous source (equation 26) cannot be 
rescaled using some weighted optical depth; rather, detailed modeling would be required. 

Three high emission measure planetary nebulae are also included in the table (IC 418, NGC 
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Fig. 2. — Hyperfine multiplet near 9.9 GHz. Three hyperfine components at 9.852, 9.876, and 
10.030 GHz appear blended due to the large natural line width, 99.8 MHz, of these transitions. 
The individual line strengths are in the ratio 1:5:2. The vertical axis is the predicted antenna 
temperature for Green Bank Telescope observations of NGC 7027, corresponding to the shell model, 
listed in Table 1 as {NGC 7027}s'. The relative populations of the 2p states are assumed negligible 
(see text), and thus the 9.9 GHz lines appear in absorption. The line profiles are completely 
dominated by natural line width and thus strong Lorentzian wings are apparent. 
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7027, and NGC 6572). These objects tend to have somewhat simpler structure than the more 
complex HII regions. Nevertheless, NGC 7027 has a well documented shell structure. In this case, 
a number of emission line diagnostics indicate rig ~ 5 x 10'^ cm~^ (Middlemass 1990). Subarcsecond 
radio images show that the emission originates from a shell with a peak emission measure of about 
2.7 X 10^ pc cm~^ (Bryce et al. 1997). Both results are consistent with an area filling factor of 
about 30%, with characteristic emission measure of about 1.2 x 10^ pc cm~^. The higher emission 
measure boosts the line-to-continuum ratio, whereas the higher density acts oppositely due to 
coUisional de-excitation. The net result in this case is a somewhat stronger estimated line with 
cont ~ —1.4 X 10 . This case is worked out in Table 1 as {NGC 7027}^, and depicted in 
Figure 2. 

Generally, compact, high emission measure objects tend to give stronger line-to-continuum 
ratios, despite their higher densities (which result in collisional de-excitation). This is because HII 
regions tend to be loosely organized along domains of constant excitation parameter U (Habing & 
Israel 1979). Using size d as a free parameter, then E oc /cP. For £■ < 3.8 x 10^ pc cm"^, the 
HII region is optically thin at 9.9 GHz and, for S = 

~._io-2_^ (30) 



Tb,cont n4 + 1.5 

where £^8 = E/{10^ pc cm~^) and n4 = n/(10^ cm~^). Thus, in the low density limit (7x4 << 1.5) 
the line-to-continuum ratio grows in direct proportion to E, and therefore to for fixed U. Of 
course, the most compact objects have densities in excess of 1.5 x 10^ cm~^, in which case the 
line-to-continuum ratio increases with d~^^'^ oc E^^^ for fixed U, since rie oc (U/d)^^'^. For emission 
measures above w 4 x 10^ pc cm~^ an HII region becomes optically thick at 9.9 GHz and the 
advantage of increased emission measure is lost. In such cases, higher densities will reduce the 
line-to-continuum ratio. 

As discussed above, the 1.1 GHz lines are considerably weaker. Under the most favorable 
conditions of high free-free optical depth and low density, we find 

AT' 

^-RKi 3.6 X 10-^ , (31) 



Tb 



cont 



assuming, as discussed above, that S = 0. These conditions are not uncommon and exceptionally 
high emission measures are not required to achieve high optical depth at 1.1 GHz. In the case of 
M 42 the optical depth at 1.1 GHz is 1.2 and the beam diluted line strength would be about 4 mK, 
versus a continuum antenna temperature of 400 K. 



6. Conclusions 



The metastable 2si/2 state of hydrogen is likely overpopulated in HII regions. Lyman-a pump- 
ing of the 2p states is expected to be negligible due to absorption of Lyman-a radiation by dust. 
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Thus, the 2si/2 — ^ 2^3/2 transitions (9.9 GHz) are predicted to appear in absorption and the 
2si/2 '^Pi/2 transitions (1.1 GHz) in stimulated emission. Because of the short hfetime of the 
final 2p states, the width of the lines is dominated completely by intrinsic line width. In effect, then, 
the power is distributed over ~ 100 MHz of line width resulting in very weak lines. In addition, 
the power is distributed over three strongly blended hyperfine lines in each multiplet. 

Searching for the 9.9 GHz lines in high emission measure HII regions offers the best prospects 
for detection. In the optically thin limit, the line strength varies as the square of the free-free optical 
depth. Predicted line-to-continuum ratios (in absorption) range up to several tenths of a percent in 
W58A, including the effects of line blending. With the Green Bank Telescope, the predicted peak 
absorption line strength may reach AT^ f» — 170 mK in this case, allowing for the redistribution of 
line strength over the three hyperfine lines. Other high emission HII regions are expected to show 
somewhat weaker 9.9 GHz lines, for example, with line-to-continuum ratios of about 0.1 percent 
and line strengths of tens of mK with the Green Bank Telescope. These predictions are uncertain, 
however, owing to biases inherent in estimating emission measures as well as selection effects in 
various inter ferometric surveys of compact HII regions. 

These conclusions apply to thermal sources. An important extension of this work would 
consider the broad line regions of active galactic nuclei and quasars in which a strong nonthermal 
microwave radiation field could influence the populations of the 2s and 2p levels, as well as provide 
a background for line absorption or stimulated emission. 

In general, detection of the fine structure lines of hydrogen will be challenging due to the 
extraordinary line width and blended structure. The observations will require meticulous baseline 
calibration and subtraction. 
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Table 1. HII Regions — Estimated line strengths. 
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0.023 


-1.0 X 


10-3 


-25. 


-1.7 X 


10-5 


1.0 X 


10-5 


0.03g 


2 


{NGC 7027}s* 


5.0 X 10* 


1.3 X 


108 




17."3 


0.34 


0.012 


-1.4 X 


10-3 


-30. 


-8.6 X 


10-5 


5.4 X 


10-5 


0.05g 


7 



=-ShelI model for NGC 7027. See text. 

Note. — Antenna temperatures (ATa) include the effects of beam dilution with the Green Bank Telescope at 9.9 GHz and 1.1 GHz (indicated by subscript G), 
and Arecibo at 1.1 GHz (subscript A). 

References. — (1) Webster & Altenhoff 1970; (2) Thomasson & Davies 1970; (3) Schraml & Mezger 1969; (4) Krassner et al. (1983); (5) van Gorkom et al. (1980); 
(6) Wynn-WiUiams (1971); (7) Bryce et al. 1997 



